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(54) Multislice magnetic resonance imaging using an array of receiving coils 



(57) In magnetic resonance imaging apparatus hav- 
ing a resistive electromagnet with a bore 6, r.f . excitation 
pulse substantially simultaneously excites a number of 
slices A to D which are phase encoded and frequency 
encoded in the usual way, the direction of the main field 
being along the bore of the electromagnet in the direc- 
tion z. The receive coil consists of an array of coils 5a 



to 5d which view different spatial regions of the imaging 
volume and the outputs of which are combined in differ- 
ent ways in order to reduce the number of slice encoding 
steps in the z-direction needed to distinguish between 
the slices A to D. Each coil of the array 5a etc. can form 
part of a two dimensional array in order to reduce the 
number of phase encoding steps in the phase encode 
direction. 
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Description • ' - . . ■ ■ ■ r ^ - 

[0001] This invention relates to magnetic resonance imaging. " • i . * 

[0002] In magnetic resonance imaging, the subject to be imaged is positioned in a strong magnetic c field, and magnetic 
s resonant (MR) active nuclei such as protons in hydrogen atoms align parallel and anffyarallel to the magnetic field, 
precessing around the direction of the field at the Larmor frequency.- - - l ^ " » ~- 

[0003] "A transmit coil applies pulses of r.f: ehergyat the Larmor frequency in a direction orthogonal* to the main field 
to excite processing nuclei to resonance, which results in the net magnetisation of all MR active nuclei being flipped 
from the direction of the main magnetic field into a direction having a transverse component in which'it can be detected 
'0 by the use of'a receive coil. " : . : *.<' - : * .... , - : f 

[0004] ^ The received signal can be spatially encoded to produce two dimensional information about the distribution 
of MR active nuclei and hence of water and' tissue. ' " : * r ' ■*=>.'■•.'■ v 

[0005] Referring to Figure' f; the main field is in the direction of the z-axis, andlhe strength of the magnetic field is 
made to vary along the direction of the z-axis" by switching on z-gradient coils. The Larmor frequency of MR active 
nuclei will then also "vary along the z -axis/arid excitation can be confined to a slice such'aV slice A^by appropriate 
choice of the bandwidth of the r.i. pulse! ' ' v ' " : J '"' ^ ♦ ; 

[0006] * Spatial information within the plane of the slice can be "obtained by reading out the'signals picked up by the 
receive coil in the presence of a gradient which varies e.g. in the y-direction, : and by turning on and turning off a gradient 
in the orthogonal direction (in this case the x-direction) before the 'readout for a number of different increments of 
20 gradient strength and polarity in order lo phase encode the received signal in that direction. — * 

[0007] There are many occasions e.g. imaging the brain of a patient, when it is desired to produce information on 
several slices A, B, C etc. in which case the procedure outlined for slice A is repeated for slice B and then for slice C etc. 
[0008] The signal picked up by the receive coil (the free induction decay - FID signal) could'appear after as little as 
a very few ms (e.g. 2 to 3ms) after the r.f. excitation r pulse, and the collection of data corresponding to one tf. pulse, 
25 including time for application of gradients and signal conditioning, coula 1 be as little as 5 to 6ms. : However, a desired 
.repetition time (TR) at which the longitudinal and transverse components of magnetisation could have recovered to 
'values necessary to produce a desired contrast could be 200ms or much longer (up to a few seconds). Therefore, if 
TR were to be chosen as 40ms, a FID signal at a first phase encode increment would be collected from slice A; and 
40ms later a second FID signal at a first phase' encode increment Would be collected from slice B. After five slices," ah 
30 FID signal of a second phase encode increment would be collected from slice A. 

[0009] In the case of a lower power magnet/a satisfactory signal-tb-noiscratiomight be produced only by completing 
the scanning procedure for all phase encode increments for all slices^ahd then repeating the data collection procedure 
a number of times. ' - * * , ' : v ' * •* : - 5 • • ■ ' > *• .-v ■ '»-..• . . / 

[0010] To iHcrease the sighal-to-noise ratio without repetitions, or with a reduced number of repetitions, a technique 
is avaitable'for exciting a number of slices e.g. four, simultahebusly/ ! usihg only one" r.f: pulse- which'behaves like four 
separate "r.f. pulses offrequehcies corresponding to the 'four slices in the presence of a magnetic field gradient, ; so that 
four times as much information can be collected after each r.f. pulse. ' 1 : '-.*-. : - r . 

[0011] Of course, each FID signal collected would now relatWto four-slices, and the slices-would have f to be distin- 
guished between in one way or another. One possible way would be to vary the r.f. phase so that e.g. on one excitation 
t - he f r ; f ; pulse. corresponding to one slice was in anti-phase compared to that corresponding to another. This set of data 
could then be identified from the signal after another excitation in which all Were in the same phase.' Permutating and 
combining such phase changes permits identification of the individual signals from the separate slices. 
[0012] Another way of distinguishing between the slices would be to apply an additional slice encoding gradient. 
After each r.f. pulse, another gradient" in the z-direciion could be turhed on and off before readout of the FID signals. 
For each phase encode gradient in the plane of the slices, a series of slice encode gradients would be applied over a 
range of gradient strengths and polarities. 

[0013] While a useful amount of information is collected after each r.f. pulse in the interests of good signal-to-noise 
ratio, the time taken for such a data collection procedure increases. 

[0014] The Applicants have appreciated that the time taken for data collection' necessitated by slice encoding can 
be reduced by the application of a technique which has been proposed for phase encoding in the image plane. 
[0015] This technique is referred to as Simultaneous Acquisition of Spatial Harmonics (SMASH): Fast Imaging with 
Radio Frequency Coil Arrays, Daniel K Sodickson, Warren J Manning, MRM 38:591 -603 (1997) and WO-A-98/21600. 
[001 6] It will be recalled that the purpose of phase encoding e^g. 'in the x -direction in Figure 1 , is to have some means 
of distinguishing the contributions to the FID signal from processing protons at various points along the x-axis. To do 
this, a gradient in the x-direction is turned on and turned off, creating a phase difference between the'precessing protons 
at different points along the x-axis, before the readout takes place. A series of gradients are in fact applied, and the 
detected signals can be processed to sort out the locations of the precessing protons which contributed to the received 
signal. (The readout gradient resolves the contributions from the protons in 'the y Section). 
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[0017] Referring to Figure 2, consider for example a coil 1 (shown dotted) for producing a saggital (vertical longitu- 
dinal) section 2 through a spine 3. In this case, the direction of the main magnetic field is along the z-axis, and the 
slices such as 2 are frequency selected by a magnetic field gradient along the x-axis. The slice 2 is spatially encoded 
by phase encoding in the z-direction and frequency encoding in the yndirectbn. A series of gradients in the z-direction 
5 are turned on and turned off before the FID signal is detected to accomplish phase encoding. At each gradient, a phase 
difference will be generated between protons precessing at the.same frequency but located at different positions along 
the length of the spine. £>nce. this is done, it is possible to calculate the, contribution to the FID signal from different 
positions along the length i of the spine. - 7 ^ '.'-■.■.;-*■■••■* . 

[0018] The SMASH technique.- refers tQan.array of. overlapping i coils 4a, 4b f 4c, 4d, £e,_4f,. 4g, 4h instead pf one 
'0 large coil 1. Each coil of the array now only viewspart of the region being imaged. The ..response patterns of the 
individual coils are shown Jn Figure 3a. — • • - ^ - • . - : -v- - v- ■> - • •■ ' 

[001 9] If the outputs of the individual coils of the array are summed, the response is ( as shown in'. Figure 3b/ and is 
similar to.that of the conventionale ingle coil. If the full amplitude output of the coil f is. summed. with the outputs of coils 
b. g, and e at reduced amplitudes, and the output of coil c and,.witb reduced amplitude, of coils b,.d, and h are subtracted, 
'5 the overall response otthepojls will be such as to create a phase, difference in the z-direction, as shown in Figure 4b. 
I hus, the FID signal wiil be different as regard to the contributions oi protons precessing in phase located at the ends 
ot the spine and the middle of the spine. The array corresponds to about one and a half.cycles when combined in this 
wny. If the array corresponded to one* whole cycle, the response would be first order, etc. 

[0020] The outputs of the coifs of the array can be combined by alternately summing and subtracting them to produce 

20 h higher order response as shown in Figure 5b. 

[0021] The SMASH technique thus partially replaces gradient phase encoding in the plane of the slice 2 by a spatial 
encoding procedure that relies on the fact that, the coijs of the ( array are sensitive to different parts, of the volume to be 
imaged. Different combinations of the outputs of those coils.produce amplitude modulation of the FID signal received 
nlong one axis oi the slice, whereas conventional phase encoding. modulates the FID signal along that axis by applying 

25 diflcrcnt magnetic field gradients to produce different phase shifts. , , 

( [0022], An advantage of the SMASH technique is that two or more output combinations can be processed in parallel 
\ox any applied phase encoding. gradient, thus reducing the number of phase encoding gradient steps required and 
speeding up the collection of data. , . _ 

[0023] Nevertheless, the SMASH technique onjy. contemplates arranging the coil array to encode Jhe MR signals in 
30 ihe image plane. . . » , . . • - • \ - - ,:. 

[0024] The invention provides magneticresonance imaging.apparatus comprising means for imparting a gradient to 
the main magnetic field. of the.apparatus, r.f. means for exciting nuclear magnetic resonance substantially sirnultane- 
ously in a plurality of slices orthogonal to the direction of the gradient in a region to be imaged, an array of receiving 
..coils for .collecting the data from, the excited slices, and means for processing the collected data using the signals from 
35 the individuai-cpils for slice encoding purposes, ; in order to produce two dimensional reaj f images of the slices. i% _ 

[0025] ..The slices may be separated'along the direction of the magnetic field gradient, or they may be contiguous, 
in which case the end result is a volumetric image. ; , . 

[0026] . ( The invention will now be described in greater detail by way of example with, reference to the. accompanying 
drawings, in which: *. , „ t . - . . ■ - - - ..- > r v.-. .• • <•* t -. • . ;- V--, 

40 ■ • - ■* - V' . ... • ' ' . . -'.^ C'V:M:/:^>ir'ni 

, . Figure .1, is a schematic representation of. known multi-slice imaging with the slice selection gradient ..along the \ z- 

axis; ^ ^ - ... . ....... >: - 0 .. : .« ' -. v -> ( - - ^ -t :.!••-.' ■ : 

Figure 2 is a schematic/epresentation of a known saggital imaging. of a spine, both with a single receive coil and 
. with an array of/eceive coils: ,. , • .. . - . t . > r - 

Figure 3a shows the arnplitude response of the individual coils of the array of Figure 2; 

Figure 3b shows the lesponse of the array when summed in a zero order combination. ( 

Figure 4a shows the amplitude response of individual coils scaled and with polarity appropriate for one combination; 

Figure 4b shows the response corresponding to Figure 4a; , v 

55 . Figure.Sa shows the amplitude response of individual coils scaled and with polarity appropriate for a higher order 

combination; , 

Figure 5b shows the response of the higher order combination; . 
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Figure 6 is' a schematic representation of tKe multi-slice imaging apparatus according'to the invention showing the 
location -of the receive coil array; - : 

- Figure 7a represents a typical 'r.f; excitation pulse: ■ • • ' r ; ■ ," v . 

Figure 7b represents the pulses effectively applied in the'presehce of a magnetic field gradient by the r.f.< pulse of 
Figure 7a, shewn for clarity "with much less -overlap than- would be used in practice":"- > < <- , 

" / • ■ ' • ^: •< . . , - ,v ,/•: ...... . .. . . . ; > 

Figure 7c shows the slice select gradients and slice ehcc<le gradients; - ' :>.- : 'us \ . ,;*v . 

Figure 7d shows the phase encode gradients^ *■■ '■•'•*-'-*-••'" ■-.■*;'".*.*:-*• - ;;< ;i . , , . ; 

Figure 7e shows th'e-frequency encode gradient;'* 1 > ^ v-.m -.i-^.' •: :., -c 
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Figure 8 shows processing electronics connected to a typical coil of the array;-^ ' - > \ - j v ^ ^ i >r i 
Figure 9 shows an alternative form of the YeceK/e coil array forthe apparatus of Figure 6; 

■ Figure 10 shows a second alternative form for the receive coil array for the apparatus of Figure 6; ; 
Figure 11 shows a third alternative form of the receive coil array for the apparatus of Figure 6; - 
Figure 1 2A is a representation of images which would result-frbm the use of the signals from respective "singlecoils; 

• Figure 12B is a representation of the image that is produced using' the signals from all four coils; 
Figure 1 3A is a profile along image 4 of Figure T2A; ; and ' * 
- Figure 13B is a profile along image 4 of Figure 12B. ■ - : ." * -■ ., : ■ -v 

[0027] ' Referring to Figure 6; the r magnetic resonance imaging'apparatus of the invention consists of an electromagnet 
of the resistive kind (not shown), havingi'a main field along the z^axis within ah interhaTbore 6, although theinvehtiori 
inapplicable to supercondoctingelectromagnets and other types of magnet as welC In additidri/the magnetic resonance 
imaging apparatus is'provided with x, yand z-gradierit coils (not shown). The z-gradient coil produces a magnetic field 
gradient ihnhe z-directibri 'which, in conjunction with the r.f. transmit coil (riot shown) is responsible for simultaneously 
energising four slices of a region bf interest of a J patient located within the- b'bre 6 of the erectromagriet, assuming the 
slices are being selected transversely i.e.- to lie in the x,y plane. The z-gradient coll produces magnetic field gradients 
in the z -direction (thefirst part bf the slice select and slice encode gradient shown in Figure 7c)." The effect bfthe single 
r.f. pulse (Figure,7a) produced by ; an r.f . transmit coil (hot shown) in c6njurictibri ; with this z-gradienHs to excite^sub- 
stantially simultaneously four slices A, B, C, D, producing the same effect as if lour r f pulses; each frequency offset 
from each other, were sim u Ifarieou sly applied. ^ - : ■ f •: / v ; .v • . - \'.:,S 

[0028] there- are actually slight timing differences between one slice and another in order to minimise* power depo- 
sition} but these are small compared with the duration of the pulses. Of -the total pulse duration, a typical variation 
between'bne slice and the next would be 10% of that total. ; - * ' v - - * ■ i - 

[0029] Note that for the purposes of illustration, these four simultaneous pulses are shbwrV spaced apart in Figure 
7b, but the'horizontal axes bf'Figure 7a to 7e otherwise represent time. These slices are orthogonal tb'the z^axis. 
[0030] The y-gradient coils are provided for spatially encoding the imaged* slices in the y-directbn, : by 'appfying the 
frequency encode gradient shown in Figure 7e when theFID signal is received in a receive coil following excitation by 
the r.f. pulse shown schematically in Figure 7a: The gradient of Figure 7e includes a rephase portion of opposite polarity 
to the read portion, although the amplitude of the two portions can vary/ Phase encoding in the x-directiori'is produced 
by means of an x-gradient coil which steps through a number of increasing values of both polarity, as shown 'schemat- 
ically in Figure 7d. This gradient is turned on and off after the r.f. excitation pulse and before the readout of the FID 
signal during the frequency encode gradient. An excitation pulse ! (Figure 7a) may be generated for each increment of 
the phase encode gradient and the FID signal read out when the frequency encode gradient is applied. This procedure 
may be repeated for each increment bf phase encode gradient. The method- is also applicable to multi-echo methods 
such as fast spin echo and gradient recalled echos such as echo planar imaging (EPI)/ 

[0031] As described thus far, the magnetic resonance imaging apparatus is 'already known, and the" value of simul- 
taneously exciting the four slices A, B, C, D, is that four times as much information can be collected for each 'excitation 
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pulse compared to the situation where only one slice at a time was excited by the r.f. pulse. A possible repetition time 
would be 200ms and a possible time for collection of the FID signal would be 40ms, so four other groups ol four slices 
could be excited before returning to the slices described. 

[0032] Previously, for each phase encode increment in the direction x, to distinguish between the slices A, B. C, D. 

s it would have been necessary to apply an increment of a separate phase encode gradient in the z<iirection, by turning 
on and turning off this .gradient,(the second lobe of the slice select and slice, encode gradient - Figure 7c), between the 
r.f. excitation pulse and readout, and this, would. have meant that the -data collection took a long time. . 
[0033] In accordance with the invention, the receive coil consists of an array of coils 5a ( 5b, 5c, 5d, the outputs of 
which are simultaneously combined in different ways at the same time-:as*slice encode gradients are applied (Figure 

10 7c). Because of these multiple signals from the coil array 5a to 5d, there are considerably less slice encode gradients 
than would hitherto have been necessary. Typical numbers of coils 5a etc. would be in the range of from 4, to 8. preferably 
8 in the case discussed here to minimise artifacts due to the sensitive volume of a coil including too much signal from 
a region from which zero sensitivity is desirable. With four coils,, the^number of. phase .increments, in. the\slice select 
direction may be reduced by a factor of four, and with eight coils, the number of phase increments in the slice select 

is direction may be reduced by a higher -factor, say, five or six. . ^ ; , : , = * - . - - ; 

[0034] Thus, data is collected in the following way. An excitation pulse is applied and the four slices A to D are 
simultaneously excited. One increment of the phase encode gradient is turned on and turned off, and the FID signal 
is read out with the frequency encode gradient applied. This procedure takes place at a first increment of the slice 
encode gradient. The procedure is repeated now at a second increment of the slice encode gradient. The procedure 

20 is repeated until all the increments of the slice encode gradient are stepped through, whereupon a second increment 
of the phase encode gradient is applied, and the slice encode gradients are again stepped through. For each data 
acquisition at an increment of the slice encode gradient, two or more separate combinations of the outputs of the array 
coils 5a to 5d are produced. These may be for the zeroth order, the Jirst order and the second order (corresponding to 
a flat response along the array, one cycle along the array, two cycles along the array, respectively). Each data acquisition 

25 therefore produces more than one line of k-space, and an array of k-space data is built up for each slice A, B, C, D. 
These are then Fourier transformed to form real two dimensional images of the four slices in the usual way. However, 
the processing of data can take advantage of any of the techniques used for. SMASH processing as described in the 
prior art referred to. 

[0035] Each coil 5a etc. is typically rectangular (Figure 8) and the coils are overlapped to reduce mutual inductance 

30 between them. Referring to Figure 8. each coil has a tuning capacitor Ct connected across it, the output of which is 
amplified in pre-amplifier 7 and mixed to an IF frequency in IF mixer 10, being bandpass filtered 8,9 before and after* 
The signal is then amplified- in amplifier 11 and digitised in analogue-tordigital converter 12^ which feeds digital signal 
processor, 13 which performs a . buffering function and in which decimation takes place. The output then passes to a. 
computer bus. A further processor (not shown) combines the outputs of the coils in the*ways described. The, amplitude 

35 and phases of-the outputs of the coils are adjusted so that the array of coils produces overall sensitivity patterns which 
correspond-to the various-spatial frequencies some of which .are: shown Jf^ Figures 3b, 4b .and 5b.; : K:Spacedata is- 
acquired, and Fourier Transformation takes. place to produce, real images of the. slice.-.. ■. t . : f - ,•■ ; . ^ 
[0036] The coils are s,uch that their sensitive volumes are desirably Gaussian or Quasi^Gaussian.^The basic restric-. 
tion on the coils .involves ensuring, that the M fields of. view of the. coils do not extend too : far,. rf sioceigeneration -of the; 

40 desirable profile patterns, becomes .more, difficult.. . . - : . ; .... y v.lu • M : ! : */ v <; iv- 

[0037] A modification is shown in Figure 10. In order to achieve the object of, ensuring tha\,fields,of view oflhe.coil.s* 
do not extend too far, the coils may be of lesser dimension along the .Z/axis with.the/esul^thatjthey then danot overlap-.; 
This would result in significant mutual inductance between the coils of the array whjch would bejundesirable. To avoid 
this, shields 16a to 16e are provided, and this reduces the-mutual inductance between the array coils 5a4o 5a> The 

45 shields are : loops tuned to the same frequency as the. coils 5a to 5d, and are shown in Figure TO at right angleslo the 
plane of the drawing. The shield loops have tuning capacitors, ^ich are tuned as a.set with the, respective, array coils.- 
The shield loops do not have to connect to any part .of the circuit of Figure 8. , ; - . 

[0038] Referring to Figure 9, it is nol necessary for there -to be a single set of coils 5a to 5d in a linear direction. A 
plurality of further sets, only one of which 14a to 1 4d is shown, can also be provided extending in the x-direction. More 

50 than two sets are required in practice. The further sets are also overlapped in the z-direction, but are also overlapped 
as well in.the x-direction, again, to,avoid mutual coupling between the coils. The outputs of coils 5a, 14a etc. are 
combined in different ways using the SMASH technique to reduce the number of phase encoding gradients, thereby 
combining the advantages of the invention with the phase encode advantage of SMASH. , 

[0039] According to a further modification, the two dimensional array of coils of Figure 9 may be wrapped around, 
55 the bore 6 (Figure 11 ). As before the coils 5a to 5d extend in the z-direction, but now each coil performs part of an 
array extended in the circumferential direction i.e. 5a, 15a, 25a, 5b, 15b, 25b, etc. The sets would extend right around 
the circumference of the former 17 on which they are mounted, but this is not. shown in the drawing for simplicity The 
advantage of this is again to dispense with some of the phase encode increments (Figure 7d) in the same way as for 
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Figure 9. In the present case, however, the coils used for reducing' the phase encode gradient steps' extend in the 
peripheral direction rather than in a straight line. The outputs of the coils 5a, 15a, 25a, etc. are combined in different 
ways to pYoduce various spatial sensitivities e.g. those shown in Figures 3b, 4b, 5b; to increase-the speed of the 
processing. The processing of the arrays 5a to 5d etc. in the z-direction remains as before. The x,y gradients would 
be combined to produce magnetic field gradients of the required orientations. ^ - - ^ 

[0040] Of course there is no need for the array to be restricted to four coils in length', for example more than four e. 
9. r eight "coils may be used. Equally, there is no restriction as to the number of coils in the peripheral direction or in the 
flat array of Figure 9. While the coils are shown as being not* overlapped' iri Figure 11, tfiey could be overlapped, or 
shielding coils could be usee!. Shielding coils cbuld be* used in^he Figure - : ' " - •* 

[0041]" There is no need for the slices excited to be^perpendicular tothe z-axisf they could be perpendicular to the x 
or y axes, or be. oblique. _ - " ; "* *■'''' ' ' ,iJi ^ ' lT . 1 : ' ' ^ ^ • : -i'- v -:" 

[0042] * there is no needier four slices to be excited simultaneously in'the , 2^^^ection. 'A , larger or smaller number of 
slices may be simultaneously exqited._ Nor is there need for the slices A to D to be spaced apart, theycould be con 1 
liguousrwith the result tha( the ! imaging apparatus produces a q^ ryz ^. h- : ** *> 

[0043] A further example's hbw'd escribed with reference to Figures i 2A.-1 2B; 1 3A, : 1 *3Bf of such magnetic resonance 
imaging, using an array of receiving coils for collecting data from a plurality of excited slices, the collected data being 
processed using the sig'nals from the individual coils, for slice encoding purposes. * ' * ■ 

EXAMPLE 

[0044] This approach relies on the use of a multi-coil array to provide "unique sensitivity information for each slice 
position. The resultant mixed slice images are untangled using matrix* algebra to solve a set of simultaneous linear 

equalions. .» 

[0045] While the example which fojlows relates to four slices, for example, as illustrated in Figures 6 and 10, the 
method is generally applicable to n coils and n slices. The method resolves simultaneously excited slices providing 
there is a minimum of the same number of coils as* slices arid each coil has a different sensitivity to any given slice. 
This sensitivity is spatially dependent and for n slices and n colls the total complex signal acquired in coill in a single 
pixel (Qt) is given by the equation: . \ - ■• ^ . 
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• ' " : ' S u x^+S^x2+ S 13 x 3 .S ln x n =C v , ■. • . . - ; . . : Equation 1 

where Sij is the complex sensitivity of coil i to slice j, xj is the spatially dependent complex signal from slice j, arid C is 
defined above, Similar expressions, for the other coils can be constructed leading to the soluble' set'of equations 1 ex- 
pressed in matrix form: / . * . "\. . \ .. . \ ; \ '* ~ ' " . ' "' ' * "* 
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Equation 2 , 



and in shorthand noiatton: 



[SHx].-=[C]- 



Equation 3 



ss 



The solution for [x] is given by: 



[*HS] .[C] 



"Equation 4 
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where [S]' 1 is the ..inverse of the complex sensitivity matrix [S]. The set of complex matrices [x] hotds the MR signal at 
each pixel for each slice. . ., . ... 

[0046] A multi-slice excitation was achieved (using for example the apparatus of Figures 6 or .10) by appropriate 
modification of a standard sine profiled pulse., such that a discrete set of slices could be excited. This work was per- 
5 formed with a tour-slice excitation and a four receive coil linear spine array. The slice positions were not commensurate 
with.the coil spacing as this is not a methodological requirement. 

[0047] The complex -sensitivity matrix was generated empirically from four acquisitions! Each acquisition excited a 
single transverse slice in a.large uniform phantom, one at each of the four positions at which the multitude excitation 
operated. These data then contain the. sensitivity of each coil to a single slice at, each position and [ this information 
10 provides thecomplex 4x4 array ([S] above) for each pixel .within this area covered by the uniform phantom. These data 
need only be acquired once for a given geometry of coils and slices. \~ 
[0048] Full simultaneous, four slice .data was .acquired with four cylindrical samples one placed at each* orthe. four 

slice positions. - , ; ^ ' - - : . - • *- - . >■-: - - - * - - 

[0049] The inverse complex sensitivity matrix joreach pixel was calculated, using (Gaussian elimination. This'and all 

is other image manipulation was performed using IDL (Research Systems, Colorado) running on a DEC Alpha worksta- 
tion. - , . f - ; ..... S; , . .._ . : ; ... r ; , x , 
[0050] All imaging was performed on a .05T Picker Asset scanner, the sensitivity matrix acquisitions and imaging 
acquisitions were all simple gradient echo sequences with 128x256 matrix, field of view 26cm,' f TR/TE 8.1/400 and 2 
sample averages, total acquisition time 102s. 

20 [0051] Uncorrected magnitude images of the test phantoms from a four slice simultaneously acquired data set can 
be seen in Figure 12A, (coils 1-4). In other words, image 1 in Figure 12A shows the image which would be produced 
if it was. built up using the output of coil 5a (Figure 6) only. Image 2 in. Figure 12A shows the image which would be 
produced if it was built up using trie output of coil 5b only. Similarly with images 3 and 4. Remembering that the object 
is.four cylindrical samples one placed at each. of the.four slice positions A-D (Figure 6), which are actually offset along 

2S the lengths.of the, images in Figure 12, it will be seen that coil 5d alone would produce an image of the cylindrical 
sample placed at the position of slice D, but would also produce faint images of the samples placed at slices B and C. 
This, illustrates the sensitivity of coil,5d to regions of the other slices A to C. The same is true of coils 5a to 5c, but this 
does not show up as markedly in Figure 12A. 

[0052] Pixel by pixel corrected magnitude images are shown in Figure 1 2B (coils 1-4). Figure 12B shows the image 

30 that is actually produced using the outputs ot all four coils 5a to 5d together with previously obtained sensitivity infor- 
mation: The artefacts which would have been produced by. the use of'a single coil have disappeared. 
[0053] The uncorrected images (Figure 12A) show the mixed data (the magnitude of array [C]) and the corrected 
images. (Figure 12B) have appropriately redistributed data .(magnitude of array [x]). Both sets of images have been 
windowed and levelled. in the same way. _ _ [ .. ' . '.."*'"" 

35 [0054] Figures 13A and 13B shows 'a single line prbfile through the magnitude image acquired by cbil 5tJ, Figure 
1 3A, before correction and Figure 1 3B, after correction. That is, Figure 1 3A shows the profile along image '4* (using the 
output of coil 5d only) in Figure 12A and Figure 13B shows the profile along slice D obtained using the outputs of alt 
four coils 5a to 5d and previously obtained sensitivity information. Jn Figure 13A, there are lesser peaks from the 
samples placed at slices B and C, but these have disappeared in Figure 13B. ( ^ 

40 [0055] The application of this technique for resolving four slice images with a four-coil array has been demonstrated. 
For n coils the principle holds for <n slices. The limitation, is that each coil must have sufficiently different complex 
sensitivity to each slice for the matrix equation 3 not to become ill-cdnditioned As an alteVnative to obtaining the 
sensitivity using phantoms-at each slice position in turn, the sensitivity c;an be obtained by mapping the whole volume 
ot coil sensitivity. Once this is done oblique slices in any orientation which satisfy the above criteria can be resolved 

45 assuming the coil geometry is fixed. * : 

[0056] Signal-to-noise in such processing is critically dependent on the noise in the array sensitive matrix [S]. Zero 
noise in [S] means the pure combination of uncorrelated noise with correlated signal, which would produce a net 
improvement in SNR. In reality as [S] is obtained empirically there will be noise and this contributes to error in the 
corrected data. 

so [0057] The principle application of the technique lies in the factor of n saving acquisition time for n slices. Contrast 
enhanced dynamic scanning, where temporal resolution is paramount, would benefit directly from such an approach. 
Spine imaging where arrays are routinely used would again benefit by shorter total acquisition times. 



55 Claims 



Magnetic resonance imaging apparatus comprising means for imparting a gradient to the main magnetic field of 
the apparatus, r.f. means for exciting nuclear magnetic resonance substantially simultaneously in a plurality of 
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slices orthogonal to the direction of the gradient in a region to be imaged, an array of receiving coils for collecting 
the data from the excited slices, and means for processing the collected data using the signals from the individual 
coils for slice encodingpurposes, in order to produce two dimensional real images of the slices. 

s 2. Magnetic resonance.imaging apparatuses claimed in Claim 1- including means for imparting magnetic field gra- 
dients in the slice encode direction, and the processing means is arranged to decode the slices using the magnetic 
field gradient steps in .the slice encode direction and the.different spatia! sensitivities of the individual coils of the 
array. ; 

io 3. Magnetic resonance imaging apparatus as claimed in Claim i or Claim 2, in which the receiving coils are overlapped 
in the direction of the gradients ■ v . ' ■ * 

4. Magnetic resonance imaging apparatus as claimed in Claim 1 or Claim 2, in which the receiving coils are separated 
from each other by shields. 

'5 : . 

5. Magnetic resonance imaging apparatus as claimed in Claim 4, in which the shields are loops tuned to the same 
frequency as the coils. ... . . . . . 

6. Magnetic resonance imaging apparatus as claimed in any one of Claims. 1 to 5, in which the slices are spaced 
20 apart from each other in the direction of the gradient. 

7. Magnetic resonance imaging apparatus as claimed in any one of Claims 1 to 5, in which the slices are contiguous 
to form a volumetric image. . , . • - • ~ - 

25 8. Magnetic resonance imaging apparatuses claimed in Claim 7, wherein each.coil of the array forms part of a further 
array in the phase encoding direction. 

9. Magnetic resonance imaging apparatus as claimed in Claim 8, in which the array is arranged in a sleeve-like 
formation encompassing the imaging volume. 

30 

10. Magnetic resonance imaging apparatus as claimed in any one of Claims 1 to 9, in which the processing means 
uses sensitivity data for the individual coils. ' ' ' * - 1 " - ----- , , ... 

11. Magnetic resonance imaging apparatus as claimed in Claim 10, in which the sensitivity data is arranged to be 
35 acquired when each of the slices is individually excited. 
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